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ABSTRACT 

We use three-dimensional simulations to study the atmospheric circulation on the first Earth- 
sized exoplanet discovered in the habitable zone of an M star. We treat Gliese 58 lg as a 
scaled-up version of Earth by considering increased values for the exoplanetary radius and 
surface gravity, while retaining terrestrial values for parameters which are unconstrained by 
current observations. We examine the long-term, global temperature and wind maps near the 
surface of the exoplanet — the climate. The specific locations for habitability on Gliese 58 lg 
depend on whether the exoplanet is tidally-locked and how fast radiative cooling occurs on 
a global scale. Independent of whether the existence of Gliese 58 lg is confirmed, our study 
highlights the use of general circulation models to quantify the atmospheric circulation on 
potentially habitable, Earth-sized exoplanets, which will be the prime targets of exoplanet 
discovery and characterization campaigns in the next decade. 
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1 INTRODUCTION 

The next frontier in extrasolar planet-hunting is the discovery and characterization of Earth-sized exoplanets — "exo-Earths". A par- 
ticularly promising route is to search for such planets around nearby M stars. M dwarf stars have several unique attributes that are driving 
exoplanet studies and astrobiology, as well as next-generation interferometry and direct imaging missions; they constitute at least 72% of 
nearby stars. As the least massive stars, they have the greatest reflex motion due to an orbiting exoplanet. Furthermore, the classical habitable 
(liquid water) zone around M dwarfs is typically located in the range ~ 0.1-0.2 AU, corresponding to orbital periods of ~ 20 to 50 days 
— well-matched to the capabilities of ground-based precision-Doppler surveys. With such short periods, hundreds of cycles of an exoplanet, 
with a mass of several times that of Earth, can be obtained within a decade, realizing factors of at least 10 in increased sensitivity for 
strictly periodic Keplerian signals and enabling Doppler reflex barycentric signals as small as 1 m s -1 to be recovered even in the presence of 
similar-amp litude stellar jitter and Poisson noise . Although these attributes have only recently become widely recognized by the astronomical 
community i Scalo et al]|2007 ; Tarter et aill2007 : Charbonneau et al]|2009h . many of the nearest M stars have been prime targets for scrutiny 
by leading precision-radial-velocity surveys for over a decade now. 



One of the most enticing and proximate exoplanetary systems being scrutinized is Gliese 581, with at least four exoplanets (Mayor et al 



2009) orbiting a nearby (6 .3 pc) M3V star. Two of the exoplanets annou nced by Mayor et al. I J2009I) are apparently "super-Earths" that 
straddle its habitable zone 1 Selsis et al. 2007 ). Recently, Vogt et al. 1 2O10h announced two more exoplanet candidates orbiting this star — 
one with a minimum mass of 3.1 (Gli ese 58 lg) and an orb ital distance of about 0.15 AU, placing it squarely within the habitable zone 
of its parent star. It is generally accepted I Lammer et al. 2O10t) that, for stellar masses below 0.6 Mq, an Earth-mass exoplanet orbiting 
anywhere in the habitable zone becomes tidally lockecQ or spin-synchronized within the first Gyr of its origin, such that it keeps one face 
permanently illuminated with the other in perpetual darkness. Such tidal locking will greatly influence the climate across the exoplanet and 
figures prominently in any discussion of its potential habitability. 



* E-mail: kheng@phys.ethz.ch (KH) 
f E-mail: vogt@ucolick.org (SSV) 

1 In this paper, we use the terms "tidal locking" and "synchronous rotation" synonymously, but this is strictly speaking only correct for an exoplanet on a 
circular orbit, which is our implicit assumption. 
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Figure 1. Mollweide projection of the temperature and velocity fields near the surface (P = 0.95-1 bar) of Gliese 58 lg as modelled by our atmospheric 
circulation simulations assuming tidal locking. Left: simulation with the fiducial value of the radiative cooling time used in the Held-Suarez benchmark for 
Earth. Right: simulation with the radiative cooling time lengthened by the ratio of the orbital period of Gliese 581g to the rotational period of Earth (a factor 
of 36.562). Colors denote temperature (in K) and the arrows represent the direction of the velocity field. The resolution adopted is T63L20 (192 X 96 X 20). 
The snapshots are taken at 1200 and 3000 Earth days, for the left and right panels respectively, after the start of the simulations. 



Independent of whether the existence of Gliese 58 lg is eventually confirmed, such disco veries motivate the study of atmospheric 
circulation on exo-Earths using three-dimensional simulations dShowman. Cho & Menoulboid) . which is the focus of the present paper. 
Our underlying philosophy is to use Gliese 581 g as a test bed and explore its atmospheric circulation — in th e absence of observational 
constraints, we as sume parameter values appropriate to the terrestrial atmosphere. Unlike in previous work (e.g. Jjoshi. Haberle & Reynolds 



1997; Joshi 2003), we study only the essential dynamics of the atmosphere, choosing not to model the radiative transfer and atmospheric 



chemistry, an approach which is commensurate with the quality of data currently available for Gliese 58 lg. We examine a suite of models 
both with and without the assumption of tidal locking. By systematically varying a set of physical parameters, we determine the major and 
minor parameters involved. We are primarily interested in the long-term, quasi-s table, large-scale circulation patterns — the climate — as 
opposed to the short-term temporal variations (the weather; IPeixoto & Oortll 1 984l) PI We describe our methods in $2] present our results in Sj3] 
and discuss their implications in Sj4] 



2 METHODOLOGY 



2.1 The GFDL-Princeton Flexible Modeling System 

We implement the spectral d ynamical core of the Flexible Modeling System (FMS) developed by the G eophysical Fluid Dynamics 
Laboratory at Princeton University l Gordon & Sternll9 82; Ande rson et alj2004 ; Heng. Menou & Phillippshoi 1). Dynamical cores are codes 
that deal with the essential dyna mics of atmospheric c irculation, treat radiative cooling in a simplified manner (via Newtonian relaxation) 
and omit atmospheric chemistry dHeld & Suareall9941) . The governin g equations solved are the primitive equations of meteorology, wher e 
the key assumption made is that of vertical hydrostatic equilibrium dSmagorinskvll 196^1 19641 ; I Washington & Parkinsoij|2005l : lvallisll2006l) . 
The boundary-layer friction between the terrestrial atmosphere and surface is treated using a simple linear prescription known as Rayleigh 
drag/friction. 

Following Held & Suarez 1 1994 ) and Heng. Menou & Phillippsl 1 2011 ), the first ti n i t = 200 Earth days of the simulations are discarded 
(unless otherwise stated); they are then run for t Iun — 1000 Earth days. Our results do not depend on the choice of tinit as long as it is 
chosen to be long enough that transient features due to initialization have been disregarded. We have verified this statement by performing 
multiple simulations with larger values of ti n ; t and witnessed no difference in our results. Choosing a larger value of t run allows zonal-mean 
flow quantities to be averaged over a longer period of time, thus producing smoother profiles as functions of latitude and vertical pressure but 
otherwise yielding no qualitative difference in the results. 

The numerical resolution adopted is T63L20 (192 x 96 x 20), which corresponds to a horizontal resolution of about 300 km. By contrast, 
the vertical pressure scale height is 



H 



kT 



60 km 



( T 



9p 



V200 Kj \2m H 14.3 ms- 



(1) 



where k is the Boltzmann constant, T is the temperature, m is the mass of the atmospheric mo lecule(s), q v is the accelerat i on due to g ravity 
and ran is the mass of the hydrogen atom. For comparison, we note that the simulations of Joshi. Haberle & Reynolds! i 1997 ) and Joshi 



IPeixoto & Oorll (l984) regard the weather and climate to be initial and boundary value problems, respectively. 
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Figure 2. Long-term, global temperature maps (in K) near the surface (P = 0.95-1 bar) of Gliese 581g. Left: with tidal locking. Right: a planetary rotation 
is one Earth day. The substellar point is located at = 180° and $ = 0°. These maps are averaged over 1000 Earth days, where the first 200 days of the 
simulation are discarded. 
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Figure 3. Same as Figure[2]but for the long-term, global, zonal wind maps. Contour levels are in units of m s 



(2003) use resolutions of T10L10 (32 x 16 x 10) and T21L22 (64 x 32 x 22), respectively. All of the simulations are started from an initial 
state of windless isothermality (uinit = 0, Ti n i t = 264 K) and executed with cons tant time steps of At = 900 s (i .e., ~ 10 5 time steps 
in total). The value of Ti n i t was again chosen following Held & Suarez 1 1994 ) and Heng. Menou & Phillippsl 1 2011 ). Our conclusions do 
not depend on this choice as long as quasi-equilibrium is attained. However, we note that in their simulations of hot Jupiter atmospheres, 
Thrastarson & Choi feOld) find the resulting structure of the atmospheric circulation to be qualitatively different depending on the assumed 
initial conditions. While it is likely that sensitivity to initial conditions scales with the depth of the atmosphere modelled, we consider such 
an exploration to be beyond the scope of our study. 

In their simplest form, the primitive equations collectively describe a frictionless heat engine, where no viscous terms exist to convert 
mechanical energy back into heat. This is operationally accomplished by the addition of a "h yperviscous" term, which also ensures numer- 
ical stability by quenching small-scale noise accumulating at the (horizontal ) grid scale (e.g., Shapiro 1970; Roe ckner & von StorcrJll98d : 



Stephensor 



Lty by quenching small-scale noise accumulating at the (horizontal) grid scale (e.g., Shapiro 19/0; Koeckner & von Storch 19sU; 
3 11994 : iMenou & Rauscheril2009l : iHeng. Menou & Phillippsl 1 201 ll) . Following iHeld & SuarezF (1994), we set the hyperviscous 



time scale to be 0. 1 Earth days in a fourth-order hyperviscosity scheme. 

Readers interested in more technical details of the model are referred to http : / /www . gf dl . noaa . gov/ fms. 



2.2 Held-Suarez forcing 



T he effects of stellar ir radiation and geometry, known as the thermal forcing, on the atmosphere are encapsulated in the forcing function 
Tf orce JHeld & Suarezlll994[) . In the classic Held-Suarez benchmark, the thermal forcing function is designed to reproduce the observed 
large-scale climate patterns on Earth, 

P 



To - ATep sin $ 



(2) 



where To = 315 K is the surface temperature at the equator (which we will scale down later), ATep = 60 K is the temperature difference 
between the equator and the poles, P represents the vertical pressure and $ denotes the latitude. The third term in equation I0 is a stabilizing 
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Figure 5. Global, long-term zonal (left) and meridional (right) wind maps for a simulation assuming tidal locking, where the radiative relaxation and Rayleigh 
friction times are set to be 36.562 times their original values. The maps are averaged over 1000 Earth days, where the first 2000 days of the simulation were 
discarded. 



term where AT Z = 10 K. Knowledge of the specific heat capacity at constant pressure c p and the ideal gas constant 1Z allow for the 
specification of k = TZ/c p . In the absence of observational constraints, we adopt terrestrial values for these quantities: c p = 1004.64 J kg - 
KT 1 , U = 287.04 J kg -1 KT 1 and k = 2/7. These values are plausible, since Cp ^ 10^ J kg 1 K 1 for a wide variety o f known gases and 
k = 2/7 has been adopted even for simulations of hot Jupiters (e.g., Menou & Ra uscher 2003; Heng. Menou & P hillipps 201 1). The surface 
pressure is initially assumed to be Pq = 1 bar, but we will explore variations in this parameter in S]3.5I We have initially not considered the 
presence of a tropopause, but this possibility is examined in $3.6\ 

Thermal forcing for a tidally-locked exoplanet can be mimicked by replacing the — sin 2 $ term in equat ion (O with cos(0 — 180°) cos $ 
(e.g. JCooper & Showmanl f2005. 20061: llVlenou & Rausche3l2009l ; lMerlis & Schneidenl20ld : lHeng, Menou & Phillippslboilh , such that the 
substellar point is located at <d — 180° and $ = 0° with G denoting the longitude. Equation (O can therefore be modified to become 



Tforce — 



Tq + AT EP cos (9 - 180° ) cos $ ■ 



A7- : hi [ £ 



cos $ 



(3) 



where Tq is now the surface temperature at the poles. Operationally, we find that Gliese 58 lg simulations with ATep = 60 K and assuming 
tidal locking do not come to quasi-equilibrium within the duration of the simulations (in the sense that the temperature map produced is 
discrepant from the forcing function), which is contrary to expectations because radiative cooling o ccurs quickly and the glob al temperature 
map should simply relax to the thermal forcing function (in the absence of initial background flow; Thrastarson & Chclboiol) . We therefore 
execute several simulations with different values of ATep to investigate this issue. We find that the global structure of the temperature 
and wind maps obtained are insensitive to the choice of ATep within the range 10-60 K, a point we will explicitly demonstrate in i)3.4l 
For operational reasons, we therefore select ATep = 10 K for the tidally-locked case because the simulations do attain quasi-equilibrium. 
Nevertheless, ATep is an unconstrained parameter of the system. 
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2.3 Physical parameters: scaling 



We next scale the value of To in equation 10 to one appropriate to Gliese 581g. Using the scaling, 



T oc £ 



1/4 -1/2 



(4) 



where £ is the stellar lu minosity and a is the distance from th e star, it follows that To = 278 K since Gliese 581 has a luminosity of 0.013 £® 
and a = 0.14601 AU jlvlavor et alj|2009l : h/ogt et al.ll2010h . Using this scaled value of To and ATep = 10 K, we get = 2 68 K. It is 

2007h . which is 



important to note that To is not the "equilibrium temperature" (i.e., blackbody equivalent) of the exoplanet I Sel sis et al. 



estimated to be T cq ~ 230 K (assuming a Bond albedo of 0.3, typical for Solar System objects; Vogt et ai] |2010l) . It is also important to 



note that typical estimates of T cq assume that energy is not transported from the permanent day to the night side (assuming tidal locking) 
of the exoplanet. On Earth, the equilibrium temperature is about 280 K, which is lower than the surface temperature of 315 K assumed in 
the Held-Suarez benchmark. The increased surface temperature is due to a combination of a non-zero surface albedo and the retention of 
(re-emitted) longwave radiation by greenhouse gases (mostly carbon dioxide and water vapour). Furthermore, the te mperature dif ference 
between the equator and the poles is determined by a combination of atmospheric dynamics and seasonal variations (e.g.. lVallisl2006h . so the 
simple scaling in equation I0 cannot be straightforwardly applied to ATep. Therefore, we retain ATep = 60 K in the case of a hypothetical 

Gliese 581 g with a rotational period of one Earth day. 

Only the minimum mass of Gliese 58 lg is currently known, but the dynamical stability analyses of Mayor et al. i 20091) and Vogt et al. 
— assuming co-planar orbits — restrict 1/ sini to have upper limits of 1.6 and 1.4, respectively. The mass of the exoplanet affects 
the assumed values of the exoplanetary radius, 

1/3 



and the surface gravity, 



M p 
M® 



M® 



-1/3 

^1 Re 



1/3 



2/3 



(5) 



(6) 



where 7?® = 6371 km and ry® = 980 cm s~ . If M p — 3.1 M®, then the scaling factor associated only with the mass is 1.46. If we instead 
have M p — 4.3 Me and 5.0 M®, then the scaling factors are 1.63 and 1.71, respectively. 

In the spirit of a scaled-up Earth, we initially assume Gliese 58 lg to have the same mass density as Earth (p P /p® = 1) and take the 
minimum mass to be the actual mass, M p = 3.1 M®, but explore the implications of raising the assumed mass to M p = 4.3 M® and 
5.0 M® in E]3.2| In £13.31 we vary the assumed value of p p //9® and explore the implications. We may already anticipate that the flow structure 
is mainly determined by the rotational rate (fl p ) and the details of the g as physics (the f rictional and radiative cooling times). 

The observed orbital period of Gliese 581g is 36.562 Earth days l Vogt et aklboid) . If the exoplanet is tidally-locked, then the angular 
rotational frequency is Q p — 1.989 x 10~ 6 s _1 . In £13.71 we explore the implications of varying the global radiative cooling time for a 
tidally-locked Gliese 58 lg — in essence, we are quantifying the relative importance of cooling versus advection. If we relax the assumption 
of tidal locking, then we can explore the effect of varying the rotational frequency of the exoplanet ( £|3,8l l. For example, if Gliese 58 lg has a 
rotational period of one Earth day, then Q p = 7.292 x 10~ 5 s _1 . 



3 RESULTS 

We begin by presenting results from a trio of baseline models, which we consider to be representative. Subsequently, we present a suite 
of simulations in which we systematically vary a set of physical parameters in order to determine which of them are significant. We conclude 
that the major parameters are the radiative cooling time r ra d and the rotational frequency 



3.1 Baseline models 

The left panel of FigureQ] shows the Mollweide projectior[f] of a snapshot from the simulation where Gliese 581g is assumed to have a 
mass of Alp = 3.1 M®, a mass density of p p — p®, a surface pressure of Po = 1 bar and is tidally-locked. The Newtonian cooling time is 
about 4 days, following the Held-Suarez benchmark for Earth (see H3.7\ . Since the rotational period of about 37 days is much longer than the 
radiative coolin g time, the structure o f the flow is sculpted by radiation rather than advection. The flow may also be modified by unresolved 
gravity waves dWatkins & Chdl2010l) . The relatively fast cooling time implies that the global temperature map relaxes approximately to the 
input thermal forcing function. 

While such visualizations are aesthetically pleasing, more insight is provided by looking at the temporally-averaged temperature and 
wind maps as functions of longitude and latitude — the long-term, quasi-stable climate. This is shown in Figure[2] where we contrast both the 
tidally-locked and non-tidally-locked cases. For the tidally-locked case, the permanent day side of the exoplanet is just within the classical 
T = 0°-100°C habitable temperature range. In the non-tidally-locked case where the rotational period is assumed to be equal to one Earth 

3 Pseudo-cylindrical projection of a globe which conserves area but not angle or shape. Also called the "homalographic projection". 
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Figure 6. Zonally- and temporally-averaged zonal wind speeds near the surface of Gliese 581g (P = 0.95-1 bar). Positive and negative values indicate super- 
and counter-rotation, respectively. 



day, the flow is dominated by advection rather than radiation with the temperatures at the equator hovering around a few degrees Celsius. 
The pair of global temperature maps in Figure[2]makes the point that conclusions on the exact locations for habitability on the surface of an 
exo-Earth depend upon whether the assumption of tidal locking is made (see also £13 - 8 1 ■ Even on the cold night side, the temper atures are 
comparable to those experienced in Antarctica where colonies of algae have been discovered and analyzed fedwards et al j[2004al lbh . All of 
these statements are made keeping in mind that temperature is a necessary but insufficient condition for habitability (see ij4j. 

Figures [3] and [4] show the global zonal and meridional wind maps, respectively. In the case of a tidally-locked Gliese 58 lg, fluid is 
transported across hemispheric scales at speeds ~ 1ms -1 , comparable to typical wind speeds on Earth. The wind patterns have a slight 
asymmetry from west to east due to the rotation of the exoplanet. If the exoplanet instead has a rotational period of one Earth day (and 
is not tidally locked), there is longitudinal homogenization of the winds with a counter-rotating jet at the equator and super-rotating jets 
at mid-latitude. The meridional wind map is now characterized by smaller structures. The slightly faster wind speeds recovered from the 
simulation with a rotational period of one Earth day are artifacts of assuming a higher value of ATep (60 K versus 10 K) — nevertheless, the 
global structure of the wind maps are robust predictions of the simulations, which we will demonstrate in £|3 .41 

To further explore the interplay between radiative cooling and advection, we execute another simulation where the radiative cooling 
(originally 4 Earth days) and Rayleigh friction (originally 1 Earth day) times are set to be 36.562 times their fiducial values — in essence, 
we are scaling by the ratio of the observed orbital period of Gliese 581g to the rotational period of Earth0 Basically, we are exploring the 
possibility that there is an unidentified cooling mechanism with a time scale that scales as the rotational period (assuming tidal locking). The 
key point is not in the exact value of this lengthening of the cooling time, but that it may well exceed the rotational period on a tidally-locked 
exoplanet. Due to the longer cooling time assumed, we now run the simulation for 3000 Earth days and discard the first 2000 days so as to 
attain quasi-equilibrium. The Mollweide snapshot of the temperature and velocity fields, as well as the long-term wind maps, are shown in 
the right panel of Figure Q] and also Figure|5] Since advection occurs somewhat faster than radiative cooling, zonal winds on the exoplanetary 
surface develop a stronger east-west asymmetry and there are hints of energy transport from the permanent day to the night side. The chevron- 
shaped featur e residing around the substellar point is reminiscent of that seen at ~ 0.1 bar in 3D atmospheric circulation simulations of hot 
Jupiters (e.g., Menou & Rauscheill2009l ; Heng. Menou & PhiliippsEoi 1). Trailing the feature are large-scale vortices spanning about a third 
of the hemisphere in size — their large sizes are a consequence of the Rossby deformation length scale being relatively larger due to the 
slower rotation of the exoplanet when tidally locked dCho et al J2003I) . In Figure[6] we show the zonally- and temporally-averaged zonal wind 
speeds obtained from our trio of baseline models. The models with tidal locking both have equatorial, super-rotating winds, where a longer 
radiative cooling time leads to faster speeds because of the increased effectiveness of advection. The non-tidally-locked model which assumes 
a rotational period of one Earth day has a counter-rotating wind at the equator; its faster speed is again an artifact of assuming a larger value 
of AT ep . Similar to how the Rossby length scale is increased because of slower rotation, the larger Rhines scale 1 Showman. Cho & Menoul 
2010) in the tidally-locked models results in broader jets across latitude. 

The trio of simulations presented here capture the essential physics of our simple model. In the sub-sections that follow, we will explore 
multiple variations on a theme, so as to ascertain the major and minor physical parameters involved in determining the atmospheric dynamics 
of an exo-Earth. In doing so, we will explicitly demonstrate some of the statements made in this sub-section. 



4 As an aside, we note that scaling down from Venus is an alternative, plausible approach since the Venusian orbital period is about 117 days. 
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Figure 7. Same as Figures|3]and|4]but for the tidally-locked case with M p = 4.3 Af® (top two panels) and M p = 5.0 JWg (bottom two panels). 



3.2 Varying the exoplanetary mass M v 

Assuming the exoplanets orbiting Gliese 581 to be in coplanar orbits, Mayor et al. i 20091) and Vogt et al. i 201fj) were able to set upper 
limits on M v of 4.3 M® and 5.0 M®, respectively, using stability analyses that considered 4 and 6 orbiting exoplanets. We thus explore the 
implications of varying M p /M® while keeping p v j 'p® = 1. As is evident from equations ((5} and {6), raising the exoplanetary mass to the 
(maximum) values implied by the stability analyses increases both the radius and the surface gravity, but in these cases by less than a factor 
of 2. As shown in Figure [7] it is therefore unsurprising that while varying M v yields minor quantitative differences in the global wind maps 
near the surface (P — 0.95-1 bar), the qualitative features are largely invariant to M p . 



3.3 Varying the exoplanetary mass density p v 

Varying the exoplanetary mass density is slightly different from varying the mass. Higher exoplanetary masses generally result in 
larger radii and higher surface gravities. However, a denser exoplanet has a smaller radius but higher surface gravity (and vice versa). To 
gain some intuition fo r the possible range of mass densities involved, we note that Kuiper Belt objects have typical densities ~ 2 g cm -3 
i Fraser&Brownlboid) . while the gas giants in our Solar System have densities ~ 1 g cm 3 . By contrast, the Earth has a density of about 
5.5 g cm -3 . It is unlikely that p v can vary by an order of magnitude or more. Thus, we explore the implications when p p /p® = 1/3 and 3 
(i.e., half an order of magnitude lower or higher than p®) and M p = 3.1 M®. 

Even with marked differences in the radii (6441 km versus 13397.9 km) and surface gravities (690 cm s -2 versus 2970 cm s -2 ) between 
the pair of simulations, the global wind maps in Figure [8] demonstrate that the mass density p v of the exoplanet is a minor parameter of the 
system. 



3.4 Varying the equator-pole temperature difference ATep 

As described earlier, the qualitative structure of the global wind maps is insensitive to the adopted value of the equator-pole temperature 
difference ATep. We now demonstrate explicitly that this is indeed the case. Figure [9] displays a set of global wind maps from simulations 
which assume tidal locking and for ATep = 20 K, 30 K and 60 K. While the qualitative structure of the wind maps are largely invariant to 



© 2011 RAS, MNRAS OOO.rTlfHI 



8 Heng & Vogt 



zonal wind (p p =p Eorth /3) 



meridional wind (p p =p Earth /3) 




50 











1 


Cr 



200 
9 (deg) 



meridional wind (p p =3p Earth ) 



-50 



,T T J- c 




i 

L 




:\\\ 


PI, 


J 





100 



200 
(deg) 



300 



Figure 8. Same as Figure|7] but for the tidally-locked case with p p = p©/3 (top two panels) and p p = 3p® (bottom two panels). 



the variation of ATep, the quantitative predictions differ — larger temperature differences between the equator and the poles lead to faster 
winds. It is likely that the equator-pole temperature difference is lower on a tidally-locked exoplanet due to the reduced effect of Coriolis 
deflection, which is the reason why we have not tested models with greater values of ATep. It is still possible that ATep > 60 K on 
exo-Earths in general, which will have strong consequences for the climate. 

Nevertheless, the fact remains that ATep is an unconstra ined parameter in our mod els. Until its value can be constrained by observations, 
quantitative predictions for the wind speeds cannot be made l Thrastarson & Choll2010h . 



3.5 Varying the surface pressure Po 

If Gliese 581g has an atmosphere, it may be thinner or thicker than that of Earth. For comparison, Venus has a surface pressure ~ 10 2 
bar due to its much thicker atmosphere, while the thin atmosphere of Mars produces Po ~ 10 -3 bar. It is therefore conceivable that the 
surface pressure Po spans many orders of magnitude in value for exo-Earths. 

To meaningfully compare simulations with different values of Po at the same vertical resolution (20 levels), we examine the zonal-mean 
zonal wind profil es as functions of the vertical pressure P and latitude $. We term such statistically-averaged flow quantities the "Held- 
Suarez statistics" l Held & Suarezfl994 : Heng. Menou & PhillippslEoi 1 ). Figure[l0]shows results from simulations with Po = 0.1, 1, 3 and 
10 bar. For the range of surface pressures explored, we witness no qualitative differences in the zonal-mean zonal wind profiles, which are 
characterized by super-rotating, equatorial winds. While we do not show them, we note that the global wind maps near the exoplanetary 
surface are very similar to those previously shown in Figures[7]and[8] 

While variations in the surface pressure appears to have little effect on the zonal-mean zonal wind profiles in the range of Po = 0.1-10 
bar, it is possible that the actual surface pressure may be much less (thin atmosphere) or greater (thick atmosphere) than the range we have 
examined. 



3.6 The presence or absence of a tropopause 

The terrestrial atmosphere closest to the surface may be divided into two components: the troposphere and the stratosphere. The strato- 
sphere is the barotropic component, meaning the surfaces of constant entropy are largely independent of latitude or longitude and increase 
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Figure 9. Same as Figures [7]and[8] but for three values of the equator-pole temperature difference: ATgp 
and AT E p = 60 K (bottom row). 



20 K (top row), AT EP = 30 K (middle row) 



with height. By contrast, the surfaces of constant entropy change with latitude and longitude in the troposphere, making it the baroclinic 
component. The height at which the troposp here transitions i nto th e stratosphere is the tropopause. 

To mimic the presence of a tropopause, iHeld & Suarea l 1994) implemented the following thermal forcing, 

Tforce = max {Tf orcc , Tstrat} , (7) 

where Tf orce is given by equation © and T s t ra t = 200 K is the (constant) stratospheric temperature. To further explore the variations on a 
theme associated with tidally-locked models, we use equation ((3} for Tf orcc and perform simulations with two different values of T s t r at- The 
first retains T str at = 200 K as in the case of Earth, while the second uses equation (0 to scale the value of T stra t down to 177 K. The point 
is not in the exact value of T str at adopted but rather the change it causes in the zonal-mean zonal wind profile. 

FigureQT|shows both the zonal-mean zonal wind profiles and the global wind maps produced by the simulations. The zonal-mean zonal 
wind profile is characterized by a super-rotating, equatorial jet situated at P w 0.2-0.4 bar (for Pq = 1 bar), where a lower stratospheric 
temperature corresponds to the jet residing higher up in the atmosphere simply because the tropopause is situated at a greater height in this 



© 2011 RAS, MNRAS 0Q0.[TliT5l 



10 Heng & Vogt 
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Figure 10. Zonal-mean zonal wind profiles as functions of the vertical pressure P and the latitude 5>, where tidal locking is assumed. Shown are simulations 
with different values of the vertical pressure: Po =0.1 bar (top left panel) 1 bar (top right panel), 3 bar (bottom left panel) and 10 bar (bottom right panel). 
Contour levels are given in m s . 

case. Interestingly, these differences do not translate into any qualitative differences for the global wind maps near the surface (P = 0.95-1 
bar) of the exoplanet, leading us to conclude that the presence of a tropopause does not have a major effect unless it reaches down to near the 
surface of the exoplanet. 



3.7 Varying the radiative cooling time T ra d 



Our dynamical core simulations account for radiative cooling via the addition of an extra term to the thermodynamic equation l IHeld & Suarez 



1994), 



(^forcc Tj 



+ Q 



where 



Q 



, T"rad,u 



1-trj 



Trad.d 



COS 4 $, 



(8) 



(9) 



the dimensionless pressure is a = P/ Po and the quantity <jj, is set to 0.7 following lHeld & SuarezH 19941) . In all of the simulations presented 
so far, we have kept r ra d, u = 4 Earth days and r ra d,d = 40 Earth days. In this subsection, we explore the implications of increasing both 
T"rad,u and T ra d,d by a constant, multiplicative factor. For convenience, we refer to the fiducial radiative cooling time simply as r ra d = 
Trad,u = 4 Earth days, where it is implied that r fa d,d/Trad,u = 10. The damping of low-level winds due to the drag or friction between the 
atmosphere and exoplanetary surface is treated by applying the Rayleigh damping coefficient, 



V- 



Rayleigh 



0. 



a < cr;,, 
a > oh, 



(10) 



to the velocity field. The fiducial value of the drag time is Tf r i C 
same constant, multiplicative factor to the Rayleigh drag time. 



1 Earth day. When we increase the radiative cooling time, we apply the 
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Figure 11. Simulations including the presence of a tropopause, where the stratospheric temperature is denoted by T 8 trat- The top, middle and bottom panels 
show the zonal-mean zonal wind profiles, zonal wind maps for P = 0.95-1 bar and meridional wind maps for P = 0.95-1 bar, respectively. The left and 
right columns are for simulations with T s trat = 200 K and 177 K, respectively (see text for details). Contour levels are given inm s . 



We have chosen the Held-Suarez benchmark as a baseline because of its simplicity. While there is no strict justification for extending 
the Newtonian cooling and Rayleigh drag schemes to exo-Earths, these are also the simplest schemes one can adopt to mimick the effects of 
cooling and drag. 

We now explore models where r ra a and rg-ic are increased by factors of a = 3, 10, 30 and 100 while retaining the assumption of 
tidal locking. Figure [T2l shows snapshots of the temperature and velocity fields from simulations that have attained quasi-equilibrium. As a 
increases and radiative cooling becomes less efficient, advection becomes more effective at transporting parts of the atmosphere heated by 
stellar irradiation from the day to the night side. Thus, the simple picture of segregated hemispheres where one hemisphere is the permanent 
day side and the other is the permanent night side starts to break down, at least from the view point of a temperature map. Some of these 
models are refutable, because when astronomical instrumentatio n becomes advance d enough to perform direct wind measurements of exo- 
Earths — such as was done for the hot Jupiter HD 209458b bv lSnellen etaD feoich . who measured the CO line in absorption with the line 
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Figure 12. Snapshots of the temperature and velocity fields near the surface (P = 0.95-1 bar), taken from simulations where the radiative cooling and 
Rayleigh friction time scales have been increased by factors of a = 3 (top left panel), 10 (top right panel), 30 (bottom left panel) and 100 (bottom right panel). 
Colors denote temperature (in K) and the arrows represent the direction of the velocity field. The resolution adopted is T63L20 (192 X 96 X 20). The snapshots 
are taken at 3000 Earth days after the start of the simulations. 

center blueshifted by about 2 km s" 1 — the predictions for the recorded spectrum will be different for each of these scenarios. For example, 
if advection dominates over radiative cooling and succeeds in the longitudinal homogenization of the temperature field, then the spectrum 
will have both blue- and redshifted contributions to the absorption lines obtained from transmission spectroscopy, thus cancelling out the 
spectral shift from line center. 

It remains to be determined, from both theory and observations, if the radiative cooling and Rayleigh friction times assumed in our 
models are consistent with the atmospheric chemistry actually found on exo-Earths. 



3.8 Varying the rotational frequency Q p 

The axisymmetric thermal forcing described in equation (O is valid when the assumed rotational period is less than the orbital period, 
e.g., in the second baseline model presented in ^3.11 with £l p = 7.292 x 10~ 5 s _1 (one Earth day). As the last variation on a theme, we 
explore models where the assumption of tidal locking is removed and the rotational period is assumed to be 1, 3, 10 and 20 Earth days 
(and where r ra d and rf r i c have been reverted to their fiducial values). Since the observed orbital period of Gliese 581g is about 37 days, we 
can use equation lO for the thermal forcing function. Such a parameter exploration of the rotational (and not orbital) frequency fl p is not 
unreasonable because tidal locking is a theoretical expectation based on the extrapolation of the properties of Solar System planets — which 
may not be representative of the exoplanets in the Universe at large — and not a direct astronomical observable. 



In this regime of axisymmetric thermal forcing, exoplanets may be characterized by the thermal Rossby number ( Mitc hell & Vallis 
20ld> . 

_ ^AT EP n a / g atep W n p (U) 

(2Q P R P ) 2 ~ V 287.04 I kg -1 K" 1 60 K / \ 7.292 x 10~ 5 s" 1 9290 km/ ' 

where 1Z is again the ideal gas constant. Since IZth < 1, a hypothetical Gliese 58 lg with a rotational period of 1 Earth day has a counter- 
rotating, equatorial wind, as shown in the right panel of Figure|3]and also Figure[6] When the rotational period is increased to 3, 10 and 20 
Earth days, the thermal Rossby number becomes 72. t h ~ 0.08, 0.9 and 4 respectively. In Figure[T3] we clearly see a transition from equatorial 
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zonal-mean zonal wind (R lh =0.009) zonal-mean zonal wind (R th = 0.08) 




-50 50 -50 50 

$ (deg) $ (deg) 

Figure 13. Zonal-mean zonal wind profiles as functions of the vertical pressure P and the latitude <E>, where the assumption of tidal locking is relaxed. Shown 
are simulations with different values of the thermal Rossby number: TZ t h ~ 0.009 (1 Earth day rotation; top left panel), 7£ t h ~ 0.08 (3 Earth day rotation; 
top right panel), T^th ~ 0.9 (10 Earth day rotation; bottom left panel) and Ti t h ~ 4 (20 Earth day rotation; bottom right panel). Contour levels are given in 
m s . 



counter-rotation to the beginnings of super-rotation (near th e surface; P ~ 0.8-1 bar) as !Z t h is increased — equatorial super-rotation is 
expected to occur when Kth > 1 JMitchell & Vallis Il2010u . However, even with a rotational period of 20 Earth days (bottom right panel 
of Figure 113b . we are already approaching the limit where the axisymmetric thermal forcing in equation (|2} starts to become invalid and 
therefore we do not perform more simulations with even longer rotational periods. Rather, our results in Figure [T3] are already sufficient to 
conclude that the rotational frequency fl p is a major parameter of the system. 

Finally, we note that equatorial super-rotation has previously been studied within the cont ext of non-axisymmetric thermal forcing 
i Suarez & Duffvlll992l ; [saravananll993l ; fKraucunas & Hartmannl2005l;lshowman & Polyanihoid) . ] t has also been previously studied using 
a more direct, non-axisymmetric momentum forcing IShell & Heldll2004ICho et alj|2008l) . 



4 DISCUSSION 

We have presented three-dimensional simulations of atmospheric circulation using the exo-Earth candidate Gliese 58 lg as a test bed. 
Our starting point is a dynamical model that is similar to the Held-Suarez benchmark for Earth, which is calibrated to reproduce terrestrial 
observations of large-scale climate patterns. Using this model as a baseline, we then present global temperature and wind maps which assume 
Gliese 581g to be tidally lo cked. The salient , qualitative conclusion gleaned from our study is that a "dance" occurs between the exoplanet 
and its host star. (See also Cho et al] 20081 .) When radiative cooling dominates over advection, the star successfully imprints its thermal 



signature on the exoplanet, which then resembles a sphere painted "half black and half white" (left panels of Figures[T]and[2} akin to the way 
one usually visualizes the permanent day and night sides of a tidally-locked exoplanet in a cartoon or schematic diagram. When advection 
occurs more efficiently than radiative cooling, the exoplanet leads the dance and smears out the thermal forcing imposed by the star evenly 
across longitude (right panel of Figure(2](. In the intermediate regime between these two cases, the exoplanet only partially succeeds in the 
longitudinal homogenization of temperature, an example of which is shown in the right panel of FigureQ] In this case, the simple picture of 
the permanent day and night sides becomes more involved. 

In the context of the classical 0°-100°C range of temperatures for habitability, our main finding is that the specific locations for 
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habitability on the surface of Gliese 58 lg — and exo-Earths in general — depend on whether the exoplanet is tidally-locked and how 
fast radiative cooling occurs on a global scale. A shortcoming of our approach is that we have neglected the effects of radiative transfer, 
which may have little effect on the large-scale climate patterns but will almost certainly alter the absolute values of the temperatures — 
for example, including the effects of CO2 may result in temperatures warmer than what we find (i.e., the greenhouse effect). The age-old 
question of whether an exopla net is inhabitable constitutes an active field and a broad range of approaches has been adopted by various 
researchers. ISelsis et al.l J2007I) used ID radiative-convective atmospheric models to as s ess the conditions on Gliese 581c and 581d, and find 
the latter to be potentially the first habitable exoplanet discovered. IWordsworth et all J20 1 lh also used a ID radiative-convective scheme, 
which included CO2 and H2O clouds, to examine the habitability of Gliese 581d. They find CO2 to be an important ingredient in maintaining 
global temperatures above the freezing point of water — increasing the surface gravity of the exoplanet reduces the CO2 column density, 
which results in global coolin g of its surface. They also find Rayleigh scattering to be less important for the redder spectra of M class 



Wordsworth etail 1I2OI0I) followed up with preliminary results from a three-dimensional climate model that focused on a pure CO2 



stars. 

atmosphere, using the early Martian atmosphere as a baseline. They raised the concern t hat CO2 is able to condense out in the colder 
regions of the exoplanet and may lead to (CO2) atmospheric collapse. iDvorak et alj 1 201o[) examine the interplay between the dynamical 
stability of systems with multiple stars and/or exoplanets and their habitable zones. [Lammer et al. 1 201fj) emphasize the coupling between 
the geophysical activity of the exoplanet and the prope rties of its parent star. Even the definition of the inner edge of the habitable zone 
depends on the amount of cloud cover on the exoplanet 1 Selsis et al]|2007h . Generally, these studies make the point that an exoplanet found 
outside of the classical habitable zone may not be uninhabitable — conversely, an exoplanet found within the zone m ay not be in habitable. 
Furthermore, a seemingly inhospitable exoplanetary surface does not exclude the possi bility of subsurface habitability ( Gold 1992). 

Although Gliese 581 itself is obser ved to exhibit weak Ca II H and K emission jMavor et al]|2009h , it is well known that M dwarfs 
in general are chromospherically active 1 Stauffer & Hartmannll 19861: IScalo et al. 2007 ) and it is likely that stellar activity will modify the 
exoplanetary atmosphere away from the terrestrial baseline JSelsis et al.ll2007h . It is also plausible that slower rotation results in weaker 
magnetospheres, which in turn weakens the ability of the exoplanet to shield itself from stellar irradiation — in the most extreme case, no 
atmosphere may exist at all. One may also study the effects of atmospheric chemistry, cloud cover, oce ans and a hydrological cycle on the 
atmospheric circulation of exo-Earths ( Jos hi. Haberle & Rey nolds 1997; Joshi 2003; Pi errehumbertl201 lh . but in the absence of observational 
constraints on the atmosphere of Gliese 581g we have chosen to omit these additional details. As such observational constraints become 
available, the models should evolve in sophistication to be commensurate with the astronomical data. When astronomical instrumentation 
becomes advanced enough to measure th e phase curves from exo-Earths, then the relative importance of advection versus radiative cooling 
may be quantified 1 Cowan & Agoll 2011 ) and hence a subset of our models may be verified or refuted. 



ACKNOWLEDGMENTS 

K.H. acknowledges support from the Zwicky Prize Fellowship and the use of the Brutus computing cluster (adroitly managed by 
Olivier Bryde et al. ) at ETH Zurich, as well as useful comments from Dick McCray and Rory Barnes on earlier versions of the manuscript. 
S.V. acknowledges support from NSF grant AST-0307493. K.H. is indebted to Kristen Menou for introducing him to this field of research 
and for many illuminating discussions. We thank the anonymous referee for useful comments which improved the quality of the manuscript. 



REFERENCES 

Anderson, J.L., et al. 2004, Journal of Climate, 17, 4641 
Charbonneau, D., et al. 2009, Nature, 462, 891 

Cho, J.Y.-K., Menou, K, Hansen, B.M.S., & Seager, S. 2003, ApJ, 587, LI 17 

Cho, J.Y.-K., Menou, K, Hansen, B.M.S., & Seager, S. 2008, ApJ, 675, 817 

Cooper, C.S., & Showman, A.P. 2005, ApJ, 629, L45 

Cooper, C.S., & Showman, A.P. 2006, ApJ, 649, 1048 

Cowan, N.B., & Agol, E. 201 1, ApJ, in press darXiv:1011.0428V l) 

Dvorak, R., et al. 2010, Astrobiology, 10, 33 

Edwards, H.G.M., Cockell, C.S., Newton, E.M., & Wynn- Williams, D.D. 2004, Journal of Raman Spectroscopy, 35, 463 
Edwards, H.G.M., de Oliveira, L.F.C., Cockell, C.S., Ellis-Evans, J.C., & Wynn- Williams, D.D. 2004, International Journal of Astrobiology, 
3, 125 

Fraser, W.C., & Brown, M.E. 2010, ApJ, 714, 1547 

Gold, T. 1992, Proceedings of the National Academy of Sciences, 89, 6045 

Gordon, C.T., & Stern, W.F. 1982, Monthly Weather Review, 1 10, 625 

Held, I.M., & Suarez, M.J. 1994, Bulletin of the American Meteorological Society, 75, 1825 

Heng, K, Menou, K., & Phillipps, P.J. 201 1, MNRAS, in press darXiv: 1010.12571 /3) 

Joshi, M.M., Haberle, R.M., & Reynolds, R.T 1997, Icarus, 129, 450 

Joshi, M. 2003, Astrobiology, 3, 415 



© 201 1 RAS, MNRAS 000. [711151 



Gliese 581 g as a scaled-up Earth 15 



Kraucunas, I., & Hartmann, D.L. 2005, Journal of Atmospheric Sciences, 62, 371 
Lammer, H., et al. 2010, Astrobiology, 10, 45 
Mayor, M„ et al. 2009, A&A, 507, 487 
Menou, K., & Rauscher, E. 2009, ApJ, 700, 887 

Merlis, T.M., & Schneider, T. 2010, Journal of Advances in Modeling Earth Systems - Discussion (JAMES-D), 2, 13-1 
Mitchell, J.L., & Vallis, G.K. 2010, Journal of Geophysical Research, in press darXiv:1008.1996 vl) 
Peixoto, J.R, & Oort, A.H. 1984, Reviews of Modern Physics, 56, 365 
Pierrehumbert, R.T. 201 1, ApJ, 726, L8 

Roeckner, E., & von Storch, H. 1980, Atmosphere-Ocean, 18, 239 
Saravanan, R. 1993, Journal of Atmospheric Sciences, 50, 1211 
Scalo, J., et al. 2007, Astrobiology, 7, 85 

Selsis, R, Kastings, J.F., Levrard, B., Paillet, J., Ribas, I., & Delfosse, X. 2007, A&A, 476, 1373 

Shapiro, R. 1970, Reviews of Geophysics and Space Physics, 8, 359 

Shell, K.M., & Held, I.M. 2004, Journal of the Atmospheric Sciences, 61, 2928 

Showman, A.P, Cho, J.Y.-K., & Menou, K. 2010, Exoplanets, ed. S. Seager, pg. 471-516 (Tucson: University of Arizona Press) 
Showman, A.P, & Polvani, L.M. 2010, Geophysical Research Letters, 37, L1881 1 
Smagorinsky, J. 1963, Monthly Weather Review, 91, 99 

Smagorinsky, J. 1964, Quarterly Journal of the Royal Meteorological Society, 90, 1 

Snellen, I.A.G, de Kok, R.J., de Mooij, E.J.W., & Albrecht, S. 2010, Nature, 465, 1049 

Stauffer, J.R., & Hartmann, L.W. 1986, ApJ, 61, 531 

Stephenson, D.B. 1994, Q.J.R. Meteorol. Soc, 120, 699 

Suarez, M.J., & Duffy, D.G. 1992, Journal of Atmospheric Sciences, 49, 1541 

Tarter, J.C., et al., 2007, Astrobiology, 7, 30 

Thrastarson, H.Th., & Cho, J.Y.-K. 2010, ApJ, 716, 144 

Vallis, G.K. 2006, Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation (New York: Cambridge Univer- 
sity Press) 

Vogt, S.S., Butler, R.P, Rivera, E.J., Haghighipour, N., Henry, G.W., & Williamson, M.H. 2010, ApJ, 723, 954 

Washington, W.M., & Parkinson, C.L. 2005, An Introduction to Three-Dimensional Climate Modeling, second edition (Sausalito: Univer- 
sity Science Books) 
Watkins, C, & Cho, J.Y.-K. 2010, ApJ, 714, 904 

Wordsworth, R.D., Forget, F, Selsis, F, Madeleine, J.-B., Millour, E., & Eymet, V. 201 1, A&A, in press jarXiv:1005.5098V 2) 
Wordsworth, R.D., Forget, F, Millour, E., Madeleine, J.-B., Eymet, V., & Selsis, F. 2010, ASP Conference Series, "Pathways Towards 
Habitable Planets", 430, 558 



© 2011 RAS, MNRAS OOO.fTllBI 



